Long-time tails and cage effect in driven granular fluids 
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We study the velocity autocorrelation function (VACF) of a driven granular fluid in the stationary 
state in 3 dimensions. As the critical volume fraction of the glass transition in the corresponding 
clastic system is approached, we observe pronounced cage effects in the VACF as well as a strong 
decrease of the diffusion constant. At moderate densities the VACF is shown to decay algebraically 
in time (t~ 3 ^ 2 ) like in a molecular fluid, as long as the driving conserves momentum locally. 

PACS numbers: 45.70.-n, 51.20.+d, 47.10.-g 



Strongly agitated granular fluids have attracted a lot 
of attention in recent years [l| . Most of the theoretical 
work which is based on microscopic dynamics has been 
done for either rather dilute or weakly inelastic sytems, 
generalising kinetic theory to gases of inelastically collid- 
ing particles. The velocity autocorrelation function Q 
as well as transport coefficients [3j have been calculated 
for the homogeneous cooling state, which has also been 
simulated for a wide range of inelasticities Q. 

Comparatively few studies have been performed on the 
stationary state of granular fluids in the moderate or high 
density regime. This is surprising, given the fact that the 
corresponding (elastic) molecular fluids have been stud- 
ied in great detail [5| and revealed several interesting fea- 
tures already in the dynamics of a single tagged particle: 

1. backs cattering as indicated by a negative velocity 
autocorrelation 

2. long-time tails due to the coupling of the tagged 
particle's density to a shear flow and 

3. a glass transition at a volume fraction (j) rs 0.58 
accompanied by a strong decrease of the diffusion 
constant as a precursor to structural arrest. 

It is our aim to understand which of these features per- 
tain to an inelastic gas and how they are destroyed by 
increasingly more dissipative collisions. This applies in 
particular to the glass transition, which has been conjec- 
tured to be related to the jamming transition in granular 
matter [6j]. 

Several experimental groups have measured the VACF 
in dense granular flow 0, d, S 03 EL • Positron emis- 
sion particle tracking allows to determine the location 
of a tracer particle and hence to measure its VACF in 
3-dimensional space. Thereby the VACF in the steady 
state of a vibro-fluidized bed has been measured [1] and 
was shown to exhibit strong backscattering effects. In 
2 dimensional vibrated layers, high speed cameras have 
been used to measure the VACF. Even though long- 
time tails seem to be beyond the experimental resolu- 
tion, these experiments give evidence for a nonexponen- 
tial decay [ll| . More recent experiments [l2j on 2 dimen- 
sional vibrated layers show strong caging effects like the 



development of a plateau in the mean square displace- 
ment. These effects may, however, be related to crys- 
tallization, which is observed in monodisperse vibrated 
layers [111, 111- 

Model — We investigate a system of monodisperse hard 
spheres of diameter a and mass m. The time evolution 
is governed by instantaneous inelastic two-particle colli- 
sions. Given the relative velocity g := vi — v 2 , the change 
of g in the direction n := (ri — r 2 )/ |(i"i — r 2 )| is 

(g-n)' = -£(g-n) (1) 

where primed quantities indicate postcollisional velocities 
and unprimed refer to precollisional ones. The coefficient 
of normal restitution e characterizes the strength of the 
dissipation. For real systems, e is a function of n and 
g [131 ]. Here we consider a simplified model with e = 
const. £ [0,1]. The elastic system is characterized by 
e = 1 and the sticky gas by e = 0. The postcollisional 
velocities of the two colliding spheres are given by = 
vi — 8 and v 2 = v 2 + S with S = ^-^(n ■ g)n. 

Due to the inelastic nature of the collisions, we have to 
feed energy into the system in order to maintain a sta- 
tionary state. This can be done either by driving through 
the boundaries, like shearing the system or vibrating its 
walls [3], or alternatively by bulk driving like on an air 
table or the experiments of ref. [l5| . Here we choose the 
simplest bulk driving fl6| and kick a given particle, say 
particle i, instantaneously at time t according to 

V<(t)=V 4 (t)+VDr&(t). (2) 

The driving amplitude i>Dr is constant and the direction 
is chosen randomly with (£j a ^ (t')) — ° T ijb"af30~{t — 

t') with the cartesian components Q , a = x,y,z dis- 
tributed according to a Gaussian with zero mean. In 
practice we implement the stochastic process by kicking 
the particles randomly with frequency /dp- 

If a single particle is kicked at a particular instant, mo- 
mentum is not conserved. Due to the random direction 
of the kicks the time average will restore the conservation 
of global momentum, but only on average. Momentum 
conservation is known to be essential for the appearance 
of long-time tails in elastic fluids. In fact the coupling 
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of the tagged particle's density to the diffusion of trans- 
verse shear is responsible for the long-time tail of the 
velocity autocorrelation. Momentum conservation also 
has important consequences for the static structure factor 
fl7l [l8| . To ensure momentum conservation not only for 
the time averaged quantity but at every instant of time, 
we introduce kicks of opposite direction for pairs of par- 
ticles [19]. However even this kind of driving conserves 
momentum only globally. To ensure momentum conser- 
vation even on small scales, we choose pairs of neighbour- 
ing particles and kick these in opposite directions. Our 
system is very close to an elastic fluid, in the sense that 
we provide the thermal energy "by hand". We regard 
this as a useful first step to investigate, which features 
of an elastic molecular fluid pertain to a driven inelastic 
granular fluid. 

The quantities we are interested in are the VACF of a 
tagged particle 

$(t) = (v i (i)-v i (0))/(v?(0)) (3) 

and its mean square displacement (MSD) 

Ar 2 (t) = ((r^)-r 4 (0)) 2 ) (4) 

Here (. . .) denotes an average over the random noise £i(t). 
Of particular interest is the diffusion coefficient, which 
is expected to decrease as we increase the volume frac- 
tion towards close packing. It can be obtained in two 
alternative ways, cither via the integral of the VACF 
3D = L dt$(i), or as the time derivative of the MSD 

d/\r 2 (t) 

6D = linif^oo — at Both definitions are equivalent in 
a stationary state. 

Method— In order to determine the MSD, the VACF 
and the diffusion coefficient, we performed event driven 
molecular dynamics simulations for several system pa- 
rameters: e — 0.7,0.8,0.9 and volume fractions 0.1 < 
(f) < 0.53725. To detect long-time tails it is very impor- 
tant to have good statistics, since the tails occur at times 
when the correlations are already quite small. Hence we 
chose a relatively small number of particles N = 10 4 
particles, but average each configuration over 1000 inde- 
pendent runs. Like Bizon et al.[20( we choose the driving 
frequency, fur, of the order of the collision frequency. 
The strength of the driving force is chosen depending on 
e such that the granular temperature attains the value 1. 
Crystallization of the system has never been observed in 
the simulation, unless we prepare the system in a crys- 
talline state initially, which was found to be stable in time 
only for e = 0.9 and the highest density (4> = 0.53725) 
investigated. 

Event driven simulations of dense systems with a con- 
stant coefficient of restitution are known to undergo an 
inelastic collapse. Several mechanisms have been sug- 
gested to avoid the inelastic collapse. Here we proceed as 
follows: We introduce a virtual hull of very small width 
for each sphere. Two approaching spheres then collide 
3 times: When the virtual hulls first touch each other, 




10"' 10° 10 1 10 2 



FIG. 1: The modulus of the VACF for <j> = 0.45 and e = 0.7, 
0.8, 0.9 (from top to bottom); s denotes the mean number of 
collisions per particle; the downward spikes in the curve for 
e — 0.9 indicate zero crossings of the VACF. 



there is no change in momentum; then the real spheres 
collide elastically when they touch; finally the inelastic 
change of momentum takes place when the virtual hulls 
touch upon receding. Thus the dissipation takes place 
only when the colliding particles are sufficiently sepa- 
rated, i.e. by the width of the hull. Each of these col- 
lisions is counted separately and we therefore divide the 
total number of collisions by three. In the simulations 
the width of the hull is taken to be 10~ 5 of the particles's 
diameter. The balance of energy input and dissipation 
requires Wp r » -^-T/m. We choose vm to achieve the 
same T for different e. It is convenient to use dimension- 
less units such that a = 2, m = 1 and T = 1. 

Results — Backscattering effects are expected to be 
strongest for high densities, when locally a cage has 
formed, enforcing reflection of the tagged particle collid- 
ing with neighbouring particles forming the cage. In Fig. 
[T]we show the modulus of the VACF for volume fraction 
<j> — 0.45 and for different inelasticities e = 0.7,0.8,0.9. 
For s = 0.9 and <\> = 0.45 we clearly observe oscillations 
in the VACF as a function of time measured in number 
of collisions: The VACF first becomes negative at t = 4, 
stays negative until t = 10 and becomes positive again 
before the signal disappears in the noise. The oscilla- 
tions are suppressed by the inelasticity. For smaller e 
the VACFs stay positive but show a dint, which is very 
pronounced for e = 0.8. We expect backscattering to dis- 
appear completely as the inelasticity is further increased 
due to two effects. First, in the sticky limit a tagged par- 
ticle is no longer reflected from its cage. Second, in order 
to achieve a stationary state of the same temperature, the 
driving force has to be increased for increasing inelastic- 
ity. Thereby the system is more strongly randomized and 
the cages are destroyed more frequently. 

For a higher volume fraction of <f> — 0.5 the VACF 
crosses the time axis for all investigated inelasticities to 



FIG. 2: The modulus of the VACF for (j> = 0.50 and different 
values for e. The symbols indicate negative values of the 
VACF. 



FIG. 4: The MSD for e = 0.7 and 4> = 0.1, 0.2, 0.35 (from top 
to bottom) as a function of s, the mean number of collisions 
per particle. 
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FIG. 3: The VACF for a system of inelasticity e = 0.7 and 
volume fractions <j> = 0.1, 0.2, 0.35 (from bottom to top). 

become positive again, as can be seen in Fig. O For 
e = 0.9, the range of negative VACF is wider than for 
<j> = 0.45, starting shortly after the first collision and 
lasting for about 10 collisions. For larger inelasticities, 
e = 0.8, 0.7, the range decreases. 

Due to strong backscattering effects at high densities 
it is difficult to investigate long-time tails at these densi- 
ties. Hence we plot in Fig. [3] the modulus of the VACF 
for intermediate volume fractions <f> = 0.1,0.2,0.35 and 
inelasticity e = 0.7. For a volume fraction of <f> — 0.35 
an algebraic tail is clearly visble and the exponent is ap- 
proximately —3/2 as in the molecular fluid. For volume 
fraction <f> = 0.2, the algebraic decay is shifted to larger 
times and for <p = 0.1 this tail seems to vanish in noise. 

The issue of long-time tails in elastic fluids has been 
traced to the coupling of the tagged particle's density to 
the diffusion of transverse shear. We expect the mode- 
coupling arguments to hold for the inelastic fluid as well 
as long as momentum is conserved locally, giving rise to 



diffusive relaxation of transverse momentum. Our results 
support this conclusion. 

The VACF has been studied previously for a freely 
cooling gas [2l], which is not in a stationary state 
so that the VACF in general depends on two time ar- 
guments, the waiting time tyy and the delay r. For 
2-dimensional systems a long-time tail of the form r _1 
was observed. Sheared granular gases were found to ex- 
hibit long-time tails in the VACF with an algebraic decay 
oc £~ 3d / 2 [23j in spatial dimension d = 2 and d = 3. 

Fig. [4] shows a double logarithmic plot of the mean 
squared displacements for a system of inelasticity e = 
0.7 and volume fractions 4> = 0.1 and 0.5. The ballistic 
regime can be clearly seen for up to one or two collision 
times. For larger times there is a crossover to the linear 
regime. Even for the dense systems no plateau is visible, 
although a hint of a developing plateau may be observed 
for the highest density of <fi — 0.53725. For the longest 
times the MSD grows linearly with time allowing us to 
extract the diffusion coefficient from the data. Fig. [4] 
reveals that the diffusion constant decreases roughly by 
a factor of 20, when the volume fraction is increased from 
= 0.1 to = 0.5. 

The simplest kinetic theory for granular gases is the 
Enskog approximation, which has been employed exten- 
sively for free cooling dynamics. It can be easily extended 
to driven systems |24| yielding an exponential decay of 
the VACF in the stationary state. This approximation 
implies for the diffusion coeffcient: 

13 1 1 Pf , . 

y/ir 8 na z g(a) ij 5 V m 

Here n denotes the number density and g(a) the pair 
correlation function at contact, which is usually approx- 
imated by the Carnahan-Starling formula [25] . This ap- 
proximation for the VACF is built on the the assumptions 
of molecular chaos (i. e. uncorrelated collisions and kicks) 
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FIG. 5: The diffusion coefficients reiative the their Enskog 
values; reference values for the elastic system from [26l ]. 

and a stationary Maxwcllian velocity distribution. Both 
assumptions are known to break down in granular fluids 
so that we expect to observe deviations from the Enskog 
theory. To quantify these, we plot in Fig. [5] the dif- 
fusion coefficient Dsim/DE relative to the Enskog value 
as a function of volume fraction together with reference 
values for an elastic system [26j. 

As in the elastic case, the dependence of Ds im /D E on 
the volume fraction is not monotonic, but the maximum 



is shifted to higher volume fractions as compared to the 
elastic case. The increase over the Enskog value for in- 
termediate densities is stronger while the decrease over 
the Enskog value at high densities is smaller as compared 
to the elastic case. Nevertheless we see a pronounced de- 
crease of the diffusion constant as the density of the glass 
transition in the elastic system is approached. 

Conclusion — We have investigated the dynamics of a 
tagged particle in a granular fluid, driven to a stationary 
state. Increasing the density we observe a strong decrease 
of the diffusion constant as the glass transition in the 
elastic system is approached. Cage effects are clearly 
visible at these high densities in the VACF, which was 
shown to oscillate as a function of time. As expected 
backscattering becomes weaker as the fluid is made more 
inelastic. We have shown that long-time tails exist in the 
inelastic fluid and are most clearly seen at intermediate 
densities. The issue of long-time tails is strongly tied 
to momentum conservation. Coupling to shear flow is 
one mechanism leading in the elastic case to an algebraic 
decay in time like t~ d ^ 2 which pertains to the inelastic 
system provided momentum is conserved locally. 
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